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ABSTRACT 

We present the results of an optical photometric monitoring program of 10 extremely 
radio loud broad absorption line quasars (RL-BALQSOs) with radio-loudness param- 
eter, R, greater than 100 and magnitude gj < 19. Over an observing run of about 
3.5-6.5 hour we found a clear detection of variability for one of our 10 radio-loud 
BALQSOs with the INOV duty cycle of 5.1 per cent, while on including the probable 
variable cases, a higher duty cycle of 35.1 per cent is found; which are very similar 
to the duty cycle of radio quiet broad absorption line quasars (RQ-BALQSOs). This 
low duty cycle of clear variability per cent in radio-loud sub-class of BALQSOs can 
be understood under the premise where BALs outflow may arise from large variety of 
viewing angles from the jet axis or perhaps being closer to the disc plane. 
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1 INTRODUCTION 

Even after four decades of intensive research the physical 
nature of quasar variability is an open question. The quasar 
optical variability has long been a tool used to set limits on 
the size of the emitting regions, to constrain emission mod- 
els, and to probe physical conditions close to the central 
black hole in active galactic nuclei (AGN). The microvari- 
ability phenomena found in radio-loud sub classes of AGN 
with i^J> 10, is quite common and is widely believed to be 
connected to the conditions in relativistic jets. Since radio- 
quiet quasars (RQQSOs) lack jets of significant power and 
extent, the microvariability seen in them may arise from pro- 
cesses on the accretio n disc itself, and thus could be used t o 
probe the discs (e.g. iGopal-Krishna. Sagar fc Wiitalll99&t ). 
Our current understanding of optical microvariability is that 
only about ~10 per cent of the radio-quiet quasars, which 
constitute ~85— 90 per cent of the quasar population, show 
microvariation, while rest of 10 per cent radio-loud (RL) 
sub-class shows about 35—40 per cent microvariability, when 
observed continuously for ~6 hour |Gopal-Krishna et al.l 
120031 : IStalin et alj|2004l : iGupta fc Joshill2005l ). Recently, the 
microvariability studies have been extended to understand 
the nature of the radio-quiet quasar sub-class with broad 
absor ption lines (BALs), the RQ-BALQSOs (|joshi et al.l 
2011), but still there has been a lack of systematic efforts 
towards the microvariability properties of their radio-loud 
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1 Radio-loudness parameter, R, is defined as the ratio of radio [5 

GHz] flux to the optical [2500 A] flux. 



counter part, i.e. radio-loud broad absorption lines, the RL- 
BALQSOs. 

The BAL quasars are characterized by the presence 
of strong absorption troughs, attributed to material flow- 
ing outwards from the nucleus with velocities of 5000 
to 50000 km s" 1 (|Green et all l200ll ). They constitute 
about 10—15 per cent of optically selected quasars (e.g . 
iReichard. Richards. Hall et alJliooa : iHewett fc Folt j|2003f ): 



which is generally interpreted as to represent the covering 
factor of an outflowing BAL wind. In addition, BALQSOs 
are pri marily believed to belong the class of radio-quiet 
QSOs l|Stocke et al.lll992l ). However with the advent of 
large comprehensive radio surveys, it has become clear that 
BALQSOs also constitute a significant fracti on of luminous 
radio quasar population (|Becker et al.|[2000l ). 

In term of microvariation of B ALQSOs, recent co mpi- 
lation of microvariability studies bv lCarinT et all i|2007l ) has 
pointed out that the BALQSOs may be an interesting class, 
as these shows about 50 per cent microvariation similar to 
the blazars, although their conclusion was base d on just 6 
RQ-B ALQSOs. This was further investigated by IJoshi et al.l 
(|2011| ) with three times larger sample size, and found that 
the true fraction of RQ-BALQSOs microvariation is just 2 
out of 19, which was found similar to the usual low mi- 
crovariability fraction of normal RQQSOs with observation 
lengths of about 4 hour. This result has rather provided sup- 
port for models where radio-quiet BALQSOs do not appear 
to be a special case of the RQQSOs in terms of their mi- 
crovariability properties. Further, it has also given support 
to the model, w here BAL o utflow appear to be closer to the 
disk plane (e.g. lElvisllioOol ). 
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In addition to probe the difference of microvariability 
properties among various class of AGN, there has been also 
considerable systematic efforts to investigate whether it de- 
pends on the radio-lou dness of the AGN or not. For instance, 
Ide Diego et al.l (|l998h have performed photometric optical 
observations of a sample of 17 core-dominated radio-loud 
quasars (CRL-QSOs) and 17 radio-quiet quasars (RQQ- 
SOs), to compare their microvariability properties. They 
found that microvariation in RQQSOs may be as frequent as 
in CRL-QSOs, wh i ch was also augmented by a similar result 
of H amfrez et al.l |2009l ). However, still there have been a 
lack of such systematic efforts to know whether or not this 
similarity also extend among the RQ-BALQSOs and RL- 
BALQSOs. This form the main motivation of our present 
work; to characterize for the first time microvariability prop- 
erties of radio-loud BALQSOs and compare the m with our 
earlie r study of the radio-quiet BALQSOs (| Joshi et al.l 
l201lh . using the same analysis and statistical method in 
both the cases. The investigation would also have an impor- 
tant implication in understanding the role of radio-loudness 
in microvariability properties and to get a clue whether the 
cause of microvariation in radio-loud and radio-quiet BAL 
quasars are similar or not. 

Paper is organized as follows: in Section 2, we describe 
our sample selection criteria, our observations and the data 
reductions while in, Section 3, we describe our analysis and 
results. Finally, discussion and conclusion are given in Sec- 
tion 4. 



2 SAMPLE SELECTION 

We select a sample of 10 radio -loud BALQSO from the com- 
pilation of IShen et al.l (|201ll ) (their online Table 1) based 
on Sloan Digital Sky su rvey, Data Release 7 (SPSS PR— 7; 
Abazaiian et al.l (|2009l )) quasar catalog (Schneide r et al.l 
2010al ). which satisfy two of our main selection criteria. 
First, the radio-loudness parameter, R, should be greater 
than 100. Here, the criterian chosen to limit radio-loudness 
parameter, R > 100, was based on our requirement that, 
sample size should be large enough for the sake of better 
statistics and the selected sources should have reasonably 
high radio-loudness, so as to help in detection of any ef- 
fect being mainly due to the difference in radio-loudness, 
between the two sub-classes i.e. the radio-loud BALQSOs 
verse s the earlier studied radio-quiet BALQSOs (|joshi et al.l 
l201lh . Second, the source need to be bright enough with mag 
gi < 19, so that even with a 1-2 m class telescope we could 
obtain a good enough signal-to-noise ratio (SNR) to detect 
fluctuations of < 0.02 mag with a reasonably good time res- 
olution of < 10 minute. 

We also limit the RL-BALQSOs to have absolute mag- 
nitudes Mi < —24.5, so that the flux contribut ion from the 
host galaxy can be assumed to be negligible l| Miller et al.l 
ll990T ). Our final sample consists of a total of 10 RL- 
BALQSO s, with 3 lobe and 7 core dominated r adio mor- 
phology (|Schneider et al.ll2010b1 ; IShen et al.ll201ll) . as listed 
in Table [1] The whole sample covers a redshift range of 
0.52 < z em «S 3.06. 



2.1 Photometric observations 

The photometric monitoring is carried out by using the 1.3- 
m Pevasthal fast optical telescope (hereafter 1.3-m DFOT) 
operated by Aryabhatta Research Institute of observational 
sciencES (ARIES), Nainital, India. PFOT is a fast beam 
(f/4) optical t elescope with poin ting accuracy better than 10 
arcsec RMS (|Sagar et al.ll201lT ). The telescope is equipped 
with Andor CCD having 2048 x 2048 pixels, with pixel size 
of 13.5 micron resulting in the field of view of 18 arcmin on 
the sky. The CCP is read out with 31 and 1000 kHz speed, 
with the respective system RMS noise of 2.5, 7 e~ and gain 
of 0.7, 2 e~/ADU (Analog to Digital Unit). The camera is 
cooled down thermoelectrically to —85 °C. We perform con- 
tinuous monitoring of each source for about 3.5—6.5 hour in 
SDSS— r passband at which our CCD system has maximum 
sensitivity. We observe each science frame for about 5—10 
minute, to achieve typical SNR grater than 25. The typical 
seeing during our observing runs was 1.5—3.8 arcsec. 



2.2 Data Reduction 

All image pre-processing steps, i.e. bias subtraction, flat- 
fielding and cosmic-ray removal are carried out using the 
standard tasks available in the data reduction software 
IRAF Q. The instrumental magnitudes of the comparison 
stars and the target source are obtained from the data 
by using Dominion Astronomical Observatory Photometryll 
(DAOPHOT II) software to perfor m the concentr ic circu- 
lar aperture photometric technique (|Stetsonlll987l ). Aper- 
ture photometry is carried out with four aperture radii, 1, 
2, 3, & 4 times of the full width at half maximum (FWHM). 
Utmost caution has been taken to deal with the seeing and 
for that, we have taken the mean FWHM of 5 fairly bright 
stars on each CCD frame in order to choose the aperture 
for the photometry of that frame. The data reduced with 
different aperture radii are found to be in good agreement. 
However, we adopted the aperture radii of ~2xFWHM for 
our final results as it has almost always provided the best 
signal-to-noise ratio. 

To derive the Differential Light Curves (DLCs) of 
quasar, we have selected two steady comparison stars from 
the same field, on the basis of their proximity in both lo- 
cation and magnitude to the quasar. The locations of the 
two best comparison stars for each RL-BALQSO is given 
in columns 3, 4 of Table [2] The typical colour differences 
g — r for our quasar-star and star-star pairs are close to unity 
(see column 7, Table (2). A detailed in vestigation quantify - 
in g the effect of colou r differences by ICarini et al.l (|l992h 
& IStalin et al.l (|2004l ) shows that the effect of colour differ- 
ences of this amount on DLCs will be negligible for a specific 
band. 

We employed a mean clip algorithm on the comparison 
star-star DLCs, so as to remove any outliers, i.e. the sharp 
rise or fall of the DLC over a single time bin which may 
arise due to improper removal of cosmic rays or some un- 
known instrumental cause, above the 3cr from the mean. The 
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Table 1. Properties of the observed radio- loud BALQSOs. 
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Ratio of the radio [5 GHz] flux to the optical [2500A] flux taken from SDSS DR7 jSchneider et alj|2010ah . 
Quasar type: l=core-dominant; 2=lobe-dominant (radio morphology classification following I Jiang et alj|2007h . 



data points corresponding to exposure resulting in such out- 
liers are also removed from the respective quasar-star DLCs 
as well. However, we should stress here that such outliers 
in our comparison star-star DLCs were usually not present 
and never exceeded two data points. Finally, the statistical 
analysis of microvariability is performed on the quasar-star 
and star-star DLCs, freed from such outliers, as shown in 

Figure cam 



3 ANALYSIS & RESULTS 

So far in the literature, the most commonly used statis- 
tics to quantify the intra-night optical variability (INQV) 
of DL Cs, is the '^-statistics" |Romero. Cellone fc CornM 
Il999l ). which is the ratio of observational scatter of quasar- 
star and star-star DLCs. However, Ide Diegd (|2010h has 
pointed out that the C is not a proper statistics, as the 
nominal confidence limit (i.e. C>2.57) for the presence of 
varia bility with 99% c onfidence, is too conservative. 

Ide Diegol (|201Cf l has shown that the other alterna- 
tives to the C-statistics are the one-way analysis of variance 
(ANOVA), x 2 — test and F— test, which are much powerful 
statistical tests to quantify the presence of microvariability. 
We noticed that for an appropriate use of ANOVA, the num- 
ber of data points in the DLC needs to be large enough so 
as to have many points in each subgroup used for the anal- 
ysis; however, this is not possible for our observations as 
we have typically only around 30-40 data points in our light 
curves. In addition, for the appropriate use of a \ 2 — test, the 
errors of individual data points need to have Gaussian dis- 
tribution and those errors should be accurately estimated. 
It has been claimed in the literature that errors returned by 
photometric reduction routines in IRAF and DAOPHOT 
are usually underestimated, often by factors of 1.3—1.75 
JCopal-Krishna et alj|2003l : ISagar et af]|2004l ; iBachev et all 
120051 ; iGoval et al.ll2012i ~ which make use of ax - test less 
desirable for such real photometric light curves. Therefore, 
to quantify the INOV in DLCs, we prefer the proper statis- 
tics that is re asonable to em ploy for differential photometry, 
the F-test (|de Diegdl201ol ) . 

Although the F-test is certainly better than the C-test, 
it should be noted that for the F-test to give a truly reli- 



able result, the error due to random noise in the quasar- 
star and star-star DLCs should be of a similar order, apart 
from any additional scatte r in the quasar-sta r DLC due to 
possible QSO variability. Ijoshi et alj (|201ll ) have pointed 
out that the standard F— test statistic does not account for 
the brightness difference in between quasar and comparison 
star as well as for the photometric error returned by the 
image reduction routines, which can result the false alarm 
detection(non-detection) of variability. For instance, if both 
comparison stars are brighter(fainter) than the monitored 
quasar, then a false alarm detection (non-detection) is possi- 
ble due to the very small(large) photon noise variance of the 
star-star DLC compared to t he quasar-star DLCs as demon- 
strated in Ijoshi et al.l (|201 ll ) . Although we have tried to se- 
lect our non-variable comparison star in proximity to both 
position and magnitude to the quasar, but it is very difficult 
to find the same for all the quasars. So here w e have used 
the "s caled F— test" statistic, as is suggested by Ijoshi et al.l 
(|201ll ), to compute our F— value as follows: 



Var(q — si) 
/i x Var(sl — s2) : 

with k, defined as, 



El 



N 2 



(q - s)/N 



Et 



Xsl-s2)/N 



Var(q - s2) 
k x Var(sl — s2) 



(° 2 (q-s)) 

(a 2 (sl- S 2)) 



(1) 



(2) 



where cr1 err (q - 



s) and a i err (sl — s2) are, respectively, 
the errors on individual data points of the quasar-star and 
star-star DLCs as returned by the DAOPHOT routine. Here, 
the term k is a scaling factor, applied to the variance of the 
star-star DLC. Firstly, this scaling factor will take care of the 
difference in magnitude between the QSO and star in quasar- 
star and star-star DLCs. Secondly, it is free from the problem 
of uncertain error underestimation by DAOPHOT/IRAF 
routines reported by many other authors, because our scal- 
ing factor depends on the ratio of the averaged squared er- 
rors and hence any factor with which errors are either under- 
estimated or overestimated will be canceled out. Therefore, 
we report our final results based on this scaled F—test (de- 
fined by F a ). 

We also note that, in standard F-test, the F-value is 
a ratio of the two sample variance, with the assumption 
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Table 2. Positions and magnitudes of the radio-loud BALQSOs and the comparison stars. 



IAU Name Date R.A.(J2000) Dec.(J2000) g r g-r 

dd.mm.yy (h m s) (° ' ") (mag) (mag) (mag) 

(1) (2) (3) (4) (5) (6) (7) 
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that error-bar in both the samples is of similar order. In 
other words, more specifically when divided by a 2 the vari- 
ance will be distributed according to \ 2 distribution, and 
hence F-distribution by definition consists of the ratio of 
two x 2 distributions, which is finally used to derive its prob- 
ability distribution. Our scaled F—test, with « term, also 
amount to the ratio of two \ 2 distribution by assigning 
a 2 = ^2o- 2 err /(N — 1) as averaged error. The only dif- 
ference here, compared to the standard F-test, is that the 
error-bar appearing in the \ 2 numerator (i.e. of quasar-star 
DLC) and denominator (i.e. of star-star DLC) do differ due 
to any brightness difference between quasar and its compar- 
ison stars. Therefore, in scaled F—test unlike the standard 
F-test, a 2 term of the numerator and denominator though 
does not exactly get canceled out, but still do preserve its 
genuine F-distribution. 

F values computed by Eq [TJ are compared individually 
with the critical F value, F^ sl/ss , where a is the signifi- 
cance level set for the test, and vqs and vss are the degrees 
of freedom of the quasar-star and star-star DLCs, respec- 
tively. The smaller the a value, the more improbable that 
the result be produced by chance. For the present study, we 



have used the significance levels, a = 0.01 and 0.05, which 
corresponds to confidence levels of greater than 99 and 95 
per cent respectively. If F is larger than the critical value, 
the null hypothesis (i.e. no variability) is discarded. Thus, 
for variability status we prefer to compare the F— value cor- 
responding to quasar-starl and quasar-star2 DLCs (i.e. F* 
& F£) from Eq. [TJ separately with the critical F value, as 
if one DLC indicates variability and other doesn't, these 
mixed signals bring into question the reality of the putative 
variability. So, a quasar is marked as Variable ('V') for a 
F- value ^ F" c (0.99), which corresponds to a confidence level 
^ 0.99; Probably variable ('PV') if the f-value is between 
F c (0.95) and F c (0.99); Non variable ('NV') if the F- value 
is less than F c (0.95). Finally, we assign a source as; of class 
variable (i.e. 'V') if both the quasar-starl and quasar-star2 
states it as 'V; of class probably variable (i.e 'PV') if either 
both DLCs show it as 'PV, or one of it as 'PV and other 
as 'V; of class non variable (i.e. 'NV') if any one of the 
quasar-starl or quasar-star2 DLC show it as non variable 
(i.e. 'NV'). 

The DLCs of our radio-loud BALQSOs sample are 
shown in Figure [TJ [2] and the corresponding results of our 
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Table 3. INOV results for the radio-loud BALQSOs. 



RL-BALQSO Date T N scaled F-test Variability?^ y^, y/i^err) 

dd.mm.yy hrs Ff,F| F c (0.95) F c (0.99) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



J004323 


42-001552 


4 


14 


11 


2012 


6 


73 


38 





76,0 


94 


1 


73 


2 


18 


NV(NV, NV ) 


3 


56 





03 


J082231 


53+231152 





28 


01 


2012 


1 


04 


33 


2 


35,2 


44 


1 


80 


2 


32 


V(V ,V ) 


1 


06 





02 


J082231 


53+231152 





11 


11 


2012 


1 


21 


23 





76,0 


93 


2 


05 


2 


78 


NV(NV, NV ) 


2 


31 





02 


J085641 


56+424253 


9 


22 


02 


2012 


3 


51 


25 


1 


40,1 


13 


1 


98 


2 


66 


NV(NV, NV ) 


1 


67 





01 


J092913 


96+375742 


9 


23 


02 


2012 


3 


88 


11 





78,1 


08 


1 


66 


2 


06 


NV(NV, NV ) 


1 


08 





02 


J095327 


95+322551 


6 


29 


01 


2012 


1 


18 


34 


1 


13,1 


26 


1 


79 


2 


29 


NV(NV, NV ) 


1 


11 





01 


J112938 


46+440325 





24 


02 


2012 


3 


10 


10 


1 


89,1 


86 


1 


70 


2 


14 


PV(PV, PV ) 


1 


75 





01 


J115944 


82+011206 


9 


28 


01 


2012 


1 


89 


11 


2 


03,1 


92 


1 


66 


2 


06 


PV(PV, PV ) 


1 


04 





01 


J121539 


66+090607 


4 


23 


02 


2012 


1 


21 


32 





82,0 


62 


1 


82 


2 


35 


NV(NV, NV ) 


1 


77 





02 


J 122848 


21-010414 


4 


22 


02 


2012 


4 


16 


48 


1 


28,0 


86 


1 


62 


1 


99 


NV(NV, NV ) 


1 


18 





01 


J160354 


14+300208 


6 


15 


05 


2012 


3 


49 


36 


2 


35,1 


91 


1 


76 


2 


23 


PV(V ,PV ) 


1 


49 





02 



a V=variable, i.e. confidence ^ 0.99; PV=probable variable, i.e. 0.95 — 0.99 confidence; NV =non-variable, i.e. confidence < 0.95. 
Variability status values based on quasar-starl and quasar-star2 pairs are separated by a comma. 

^ Here k = (c 2 (q — s)) / (a 2 (si — s2)) (as in Eq.[2]l, is used to scale the variance of star-star DLCs for the scaled F—test. 



analysis are summarized in Table|3] In the first four columns, 
we list the object name, date of observation, the duration 
of our observation, and the number of data points (N p0 i„t s ) 
in the DLC. The column five gives the scaled F—test values 
corresponding to both quasar-starl and quasar-star2 DLCs. 
The columns six, seven mention the critical F-value for the 
corresponding 95 per cent and 99 per cent significance level. 
In column eight, we give the variability status of our source, 
based on the scaled F—test for both the quasar-starl and 
quasar-star2 DLCs. Column nine lists the square root of 
scaling factor, t/k, where k = {a 2 (q — s)/a 2 (sl — s2)} (as in 
equation [T]), and has been used to scale the variance of the 
star— star DLCs while computing the F value in the scaled 
F—test. The last column gives our photometric accuracy, 
\Z{ a i,err) m the quasar— star DLCs, which is typically of 
0.01-0.03 mag. 



3.1 The INOV duty cycle (DC) 

The duty cycle of INOV was computed as is given by 
iRomero. Cellone fc Cornbl (| 19991 ), 



y n t Ni(i/Ati) 

DC = lOO ^ 1 / percent 



(3) 



where Ati = Afi i0 (, s (l + z)^ 1 is the duration of the 
monitoring session of a source on the i th night, corrected for 
its cosmological redshift, z. Since the observing run time for 
a given source is not same for different nights observation, 
so the computation of DC has been weighted by the actual 
monitoring duration Ati on the i th night. Here, Ni was set 
equal to 1 if INOV was detected, otherwise Ni = 0. 

In 8 nights monitoring (i.e. for ~46 hour) of our 10 
RL-BALQSOs, we find a clear signature of variability for 
one RL-BALQSO resulting in INOV duty cycle (DC) of 5.1 
per cent. While on including the three sources assigned as 



probable variable ('PV, class), the DC increases to 35.1 per 
cent. 

Here, it is also worth comparing the DC of RL- 
BALQSOs with our earlier study of 19 RQ- BALQSOs for 
microvariability properties (|Joshi et al.ll20lil ). For that, we 
have computed the INOV DC value of RQ-BALQSOs sam- 
ple, by using Eq. O An INOV DC of 6.41 and 34.28 per 
cent was found respectively for clear detection of variability 
(two out of 19 DLCs) and probably variable cases (six out 
of 19 DLCs), which is almost similar to our above DC of 
RL-BALQSOs. 



4 DISCUSSION & CONCLUSIONS 

Soon after the discovery of quasars, their behavior of chang- 
ing flux in minutes to hours and amplitude from few tenth 
to a few magnitude i.e. the microvariability, was reported 
(|Matthews fc Sandagd Il963l ). Thereafter, plenty of efforts 
have been made to correlate it with several other prop- 
erties, such as quasar orientation, radio loudness and po- 
larization, driving the microvariation in various sub-classes 
of AGN. The emerging picture is that the microvariability 
does not occur as o ften in radio-quiet objects as it does in 
radio - loud sources dGopal-Krishna et a l. 2003; Stalin et alj 



12004 [c arini et al. 1 120071 ). However, a multiband monitor- 
ing study of 22 radio- quiet and 22 core do minated radio- 
loud (CRL) quasars by I Ramirez et al.l (|2009l ) . has confirmed 
that the frequency of microvariability in both these quasar 
sub-classes are same: similar to t he earlier re s ults f ound by 
Ide Diego et~all (|l99Sl ). Recently, iGoval et all (|2012T l has in- 
vestigated the role of optical polarization along with the 
relativistic beaming in microvariability property for a large 
sample of 21 core dominated quasars (CDQs), consisting of 
12 low polarization (LP) and 9 high polarization (HP) QSOs. 
With the INOV duty cycle of 28 per cent and 68 per cent re- 
spectively for LPCDQs and HPCDQs, they concluded that 
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the relativistic beaming is normally not a sufficient condi- 
tion, but high optical polarization is also a necessary condi- 
tion for strong INOV. 

In general, the microvariability seen in radio-loud 
objects is now widely believed to come from turbulence 
in a relativistic jet pointed at or nearly along our line 
of sight. However, in the case of radio-quiet objects, the 
jets are weaker than radio-loud jets or quenched near 
their origin point due t o the effect s of b lack hole spin 
l|Wilson fc Colbertl Il995l iBlandfordl l200fj| ) or magnetic 
configurations ( Meier] 2003 ). To know which mechanism 
among these various possibilities is at play for the origin 
of microvariability on different sub-classes of AGN, the 
constraints obtained by extending the microvariability 
studies to other remaining AGN sub-cla ss such as BALQ - 
SOs have given important clues (see e.g. IJoshi et aljfeoilf ). 
For instance, the early evidence from spectro-polarimetry 
of BALQSOs give an idea that they a re viewed closer to 
the d isk plane i.e. nearly edge-on (e.g. iGoodrich fc Millerl 
Il995h . It also supports the models wh ere outflow comes 
closer to the disk plane (e.g. lElvis! [20001 ) . It can be a clue 
to understand the radio-loud/radio-quiet dichotomy as the 
BAL quasars appeared to be exclusively r adio-quiet, where 
outflo ws are not well collimated. Later iBrotherton et all 
(|2006h have shown that BALQSOs are viewed at a large 
variety of viewing angles (~ 15° from the jet axis to nearly 
edge-on) . Constrained by radio flux density variability to lie 
within 35° of a relativistic jet, give rise to a new population 
of polar BALQSO i.e. those viewed perpendicular to 
the disk (e.g. IZhou. Wang. Wang. Wang. Yuan fc Lul 
20061; 



iGhosh fc Punslvl 



Montenegro-Montes. Mack. Vigotti et all 
Dpi. Kawaguchi. Kono et al.l I2009T ) . where outflows are 



2007 
200S 



confined to well-collimated relativistic jets. 

In these contexts, here we have presented our results 
of intranight optical variability observations of 10 extremely 
radio-loud (i.e. R> 100) BALQSOs, observed mostly in the 
first half of year 2012. These are the first extensive obser- 
vations of radio-loud sub-class of BALQSOs with an aim to 
probe their microvariability nature. In other words, we have 
investigated the effect of presence of BAL as well as radio- 
loudness on microvariability properties. The application o f 
proper statistics, i.e. the "scaled F—test" ll Joshi et al.ll201l] ) 
has allowed us to find the genuine confirmed INOV DC of 5.1 
per cent (one variable DLC out of total 10), which becomes 
35.1 per cent if we also include the three probable variable 
(i.e. 'PV') cases as well. This DC is free from any spurious 
cause, as could be present at times while using the standard 
F-test i.e. introduced by the brightness difference in quasar 
and its comparison star (e.g. Sectio n EH), since we hav e used 
scaled F—test, described in detail in IJoshi et alj l|201ll ) . Hav- 
ing used the sa me statistic a l and analysis method, both in 
this work and IJoshi et alj (|201uY where they studied mi- 
crovariation properties of 19 radio-quiet BALQSOs, we car- 
ried out the comparison of these two results and found that 
our resulting INOV DC of 5.1 per cent matches with their 
result of 6.41 per cent fraction (where they found out of 19 
only 2 DLCs as variable, also see Section 13. ip . This sug- 
gest that radio-loudness has not any significant role at least 
in the case of microvariability properties of BALQSOs, but 
larger sample size of radio-loud BALQSOs will help to say 
about it more firmly. 



As a general implication, if we would have also got high 
D C of microvariation in our RL-BALQSOs, such as found 
by I Zhou et ai] (|2006h & IGhosh fc Punslvl (|2007l 1 based on 
their radio monitoring, then it would have suggested that 
the BAL outflow are aligned very close to the jet axis (about 
< 35 deg) which may be responsible for the higher DC due 
to Doppler boosting. On the other hand, the common wind 
based model of BAL outflow demand the outflow to be pref- 
erentially eqjiatorial when the accretion disk is almost edge 
on (e.g. lElvisllioobl ). predicting the low DC as is also sug- 
gested by our this investigation. Though only with such 
microvariability studies, it may be difficult to distinguish 
clearly among many possibilities in a quantitative manner; 
but with our this study we certainly could conclude that the 
microvariation of both radio-loud and radio-quiet BALQSOs 
are very similar, and the level of microvariation occurrence 
in both cases is quite low as is usually found for radio-quiet 
QSOs, the theoretical implication of which for AGN unifi- 
cation scheme needs further exploration. 
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APPENDIX A: BRIEF NOTES ON 
INDIVIDUAL SOURCES 

Al J004323. 42 001552.5 

J004323.42-001552.5 is a high redshift (i.e. z=2.797) 
BALQSO with a balnk ity indej£|(BI) of 591.10 
km s _1 (|Gibson et al.l 120091 '). It appears as non- variable 
according to our scaled F—test with an observation length 
of 6.73 hr. 



A2 J082231. 53+231152.0 

J082231. 53+231152.0 is a low redshift (i.e. z=0.653), 
Mg 11 low ionizatio n broad absorption line (loBAL) quasar 
(jZhang et al.ll2010l ). This source has been monitored twice, 
firstly on the night of 28 January 2012 for about 4.04 hr 
and secondly on the night of 14 November 2012 for about 
4.21 hr. On applying the scaled F—test statistics, this source 



The Balnicity Index metric is defined as 
r.3000 r 



BI : 



■_>MMI[I 



1 



f(v) 
0.9 



Cdv. 



(Al) 



Here, the limits of the integral are in units of km s , f(v) is the 
normalized flux as a function of velocity displacement from line 
center. The constant C = everywhere, unless the normalized 
flux has satisfied f(v) < 0.9 continuously for at least 2000 km s _1 , 
at which point it is switched to C = 1 until f(v) > 0.9 again. 
Based on this definition, objects are classified as BALQSOs if 
their BI > km s _1 . 
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Figure 1. Differential light curves (DLCs) for the radio- loud BALQSOs in our sample. The name of the quasar and the date and duration 
of the observation are given at the top of each night's data. The upper panel gives the comparison star-star DLC and the subsequent 
lower panels give the quasar-star DLCs as defined in the labels on the right side. Any likely outliers (at > 3<r) in the star-star DLCs are 
marked with crosses, and the data corresponding to such flagged exposure are also removed from quasar-star DLC, for the final analysis. 
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Figure 2. Same as in Figure 1. 



has indicated a clear signature of micro-variability in its first 
DLC lasting ~4.04 hr while appeared as non-variable in its 
second DLC (e.g see Table [3]). This makes it a promising 
candidate for future INOV monitoring. 



A3 J085641. 56+424253. 9 

J085641. 56+424253. 9 is a C IV high ionization broad ab- 
sorption line ( HiBAL) quasar, with a balnicity in dex of 
820.60 km s" 1 (|Trump et alj|2006l ; [Gibson et alj|2009f ). This 
BALQSO is also identifie d as a high redshift, z=3.061, weak 
emission line quasar by iDiamond-Stanic et alj (|2009T ). We 
observed this source for ~3.51 hr duration for which the sta- 
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tistical analysis of its DLC does not give any clue of rapid 
variability. 



A4 J092913. 96+375742. 9 

This source has a massive black hole of mass 8.9xl0 9 M ( 



q, 



estim ated using its Mg II line width (see e.g. IShen et alj 
I2008T ). This source is a HiBAL with BI value corresponding 
to its C IV line is 1204.7 km s~ . We found this source to be 
non-variable during our 3.9 hr observation. 



A5 J095327.95+322551.6 

This source is a flat spectrum r adio quasar type with a 



4.8 GHz flux density of 184mJy ( Healev et al 



20071 ). This 



source is detected as BAL in iGibson et all (120091 1 SDSS 
DR-5 BAL quasar catalog, with a C IV balnicity index of 
267.0 km s _1 . This source has not shown any signature of 
microvariability, for an observing run of 4.18 hr. 



A6 J112938.46+440325.0 

This BAL quasar has 



C iv BI value of 806.5 
km s _1 (|Gibson et alj|2009h and has a massive black hole 
of, 1.32xl0 10 M(7) as estimated by using the Mg II line width 
|Shen et al.ll2008T ). This source was found to be probably 
variable during the course of our 3.40 hr observation, which 
makes it a potential source for further microvariability study. 



A7 J115944.82+011206.9 

This source, also known as, PKS 1157+014, has 
long been k nown to possess a 21 cm absorption at 
Zabs = 1-94 (IWolfe et alj|l986h associate d with the DLA 



jColbert fc Malkanl 120021+ litocke et all (| 19841 ) found it 



as unresolved point source in 488MHz VLA observation, 
with flux value of 108mJy. It is a HiBAL with balnicity 
index of 299.3, 632.9 km s _1 , respectively for Si IV and 
Civ (|Gibson et alj|2009h . This source is a promising candi- 
date for future investigation as it appears probably variable 
for our observing run of 4.89 hr. 



A8 J121539. 66+090607.4 

This source is a member of large bright quasar survey and 
is also known as LBQS 1213+0922. It has a massive black 
hole of m ass 5.12 x 10 9 Mq e stimated by using the C IV line 
width by IShen et all (|200ST >. This is a narrow trough high 
ionization BAL (i.e. nHiBAL) QSO with BI value 116.30 
km s _1 , estimated by C IV line. During the monitoring of 
this source, few exposures got corrupted (about maximum 
8-9) due to some technical problem, which we have removed 
manually resulting in few gaps in its DLC. Over the mon- 
itoring duration of ~4.21 hr, we found this source to be 
non- variable. 



A9 J122848. 21-010414.4 

This is a HiBAL quasar, with BI value of 17.0 km s~ 1 , cor- 
responding to the C IV line. It appears as non-variable ac- 
cording to our scaled F—test with an observation length of 
4.16 hr. 



A10 J160354.14+300208.6 

This is a HiBAL quasar with a massive black h ole of mass 
1.46 x 1O 9 M estimated by using the Mg 11 line |Shen et al.l 
120081 ). Its BI value is 84.7, 480.0 km s _1 , respectively for 
Si IV and C IV line. Over an observation length of 3.49 hr, 
we found it as probably variable (i.e. 'PV'), which makes it 
a potential candidate for further microvariability study. 
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